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THE  EXISTENCE  OF  GENERALIZED  EIGENFUNCTIONS 
IN  UNDERWATER  ACOUSTICS 


INTRODUCTION 

The  acoustic  field  due  to  a  point  source  in  a  horizontally  stratified 
waveguide  can  be  expressed  as  a  contour  integral.^  This  contour  follows  a 
path  in  the  complex  plane  that  separates  the  singularities  of  the  depth-  and 
range-separated  Green's  functions.  Because  it  is  assumed  that  the  depth- 
separated  problem  is  finite,  the  depth-separated  Green's  function  will  have 
poles,  but  no  branch-cut  singularities.  The  poles  of  the  Green's  function  are 
the  zeros  of  the  characteristic  equation  and  are  the  same  as  the  eigen¬ 
values.  When  the  contour  is  enclosed  around  the  poles,  the  solution  can  be 
expressed  as  the  sum  of  the  residues  at  the  poles.  If  the  depth-separated 
problem  is  self-adjoint,  then  the  poles  of  the  Green's  function  are  simple, 
the  zeros  of  the  characteristic  function  have  a  multiplicity  equal  to  one, 
and  the  residue  series  is  the  usual  normal  mode  solution.^ 

If  the  depth-separated  problem  is  non-self-adjoint,  the  poles  of  the 
Green's  function  may  not  be  simple,  the  zeros  of  the  characteristic  function 
can  have  a  multiplicity  greater  than  one,  and  the  residue  series  assumes  a 
more  complicated  form.^  The  eigenfunctions  must  be  augmented  w'.n  the 
generalized  eigenfunctions  (called  associated  functions  in  Naimark^)  to 
yield  an  expansion  theorem.  An  example  of  this  situation,  with  separated 
boundary  conditions,  is  shown  in  a  useful  book  by  Friedrian.^  His  example 
corresponds  to  an  isovelocity  waveguide  bounded  by  a  free  surface  and  a 
special  complex  impedance  bottom  boundary  condition.  Similar  examples 
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occur  in  acoustic  ducts  with  admittance  boundary  conditions.^ 

Non-self-adjoint  problems  occur  in  underwater  acoustics  when  the 
square  of  the  wavenumber  is  given  a  complex  value  to  allow  for  volume 
attenuation.  This  report  shows  that  complex  wavenumbers  can  give  rise  to 
multiple  eigenvalues  and  generalized  eigenfunctions  in  the  same  way  as 
complex  impedance  and  admittance  boundary  conditions  give  rise  to 
multiple  eigenvalues. 
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THEORY 


Consider  a  waveguide  where  the  square  of  the  wavenumber  k2(x)  = 
g  X  +  h  is  a  complex-valued  linear  function  of  the  depth  x  for  0  ^  x  <  H.  The 
waveguide  is  assumed  to  be  bounded  by  a  free  surface  and  bottom.  The 
homogeneous  depth-separated  problem  is  given  by  the  complex  Sturm- 
Liouville  boundary  value  (eigenvalue)  problem 


+[gx  +  h-  X,]u  =  0,  (1a) 

dx2 

u(0).0.  (1b) 

u(H)«0.  (1c) 


where  X  is  the  complex  separation  parameter.  A  nonzero  function  <1)^  (x)  that 

satisfies  equation  (1)  is  an  eigenfunction  and  the  corresponding  Is  an 

eigenvalue.  The  Green's  function  G(x,s,  X)  of  the  depth-separated  problem 
satisfies  the  nonhomogeneous  counterpart  of  equation  (la)  with  a  delta 
function  source  at  x  »  s.  The  Green's  function  also  satisfies  the  boundary 
conditions  in  equations  (1b)  and  (1c). 

Suppose  that  the  differential  operator  L  is  defined  by 

L  u  .  diii  +  [g  X  +  h  ]  u  (2) 

dx^ 

for  functions  satisfying  the  boundary  conditions  in  equations  (1b)  and  (1c). 

The  eigenfunction  corresponding  to  X,p  satisfies  (L-X,n)<l)p  =  0.  A 
generalized  eigenfunction  (of  rank  2)  corresponding  to  Xp  is  a  nonzero 
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function  that  satisfies  (L-Xn)2yp  .  o.  but  (L-Xn)\Kn  0.  The  Green’s 
function  satisfies  (L-X)G  «  5(x-s). 

Equation  (1a)  can  be  solved  with  Airy  functions  after  the  linear  change 
of  variables  is  made: 


z(x) 

-  -g'2/3  [  g  X  +  (h  -  X)  ]  . 

(3) 

Then  equation  (1)  becomes 

i&lL  .  z  w  ~  0  , 

(4a) 

w(zi)  -  0  , 

(4b) 

W(Z2)  -  0  , 

(4c) 

where  w(z)  -  u(x),  z-i  -  z(0),  and  Z2  »  z{H).  Grosjean  and  De  Meyer®  have 
shown  that  there  exist  distinct  complex  numbers  Zf  and  Z2  and  a  corre¬ 
sponding  nonzero  solution  of  equation  (4)  that,  in  addition,  satisfies 

f^2 

w2(z)  dz  -  0  ,  (5) 

•'zi 


where  w  is  determined  uniquely  up  to  an  arbitrary  nonzero  multiplicative 
factor. 

The  solution  of  Grosjean  and  De  Meyer®  Is  transformed  into  a  solution 
of  equation  (1)  with 
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g 


(6a) 


h3 


and 


n-2/3  . 

^  J>  ’ 

(Zl  -  22)^ 


(6b) 


h2 


(6c) 


The  resulting  solution  <j)p(x)  ■  w[z(x)l  is  an  eigenfunction  of  equation  (1) 
and,  in  addition,  satisfies 


4^(x)  dx  * 


0  . 


(7) 


Once  h  is  chosen,  the  A,  =  Xp  in  equation  (6c)  is  the  eigenvalue  correspond¬ 
ing  to  <|)n. 


If  were  a  simple  pole  of  the  Green's  function,  the  residue  at 
would  be^ 


1-1 


<|)n(x)  dx 


<l>n(x)  <{»n(s)  , 


(8) 


which  is  impossible  in  view  of  equation  (7).  Thus,  is  not  a  simple  pole  of 
the  Green's  function,  and  a  generalized  eigenfunction  must  exist  that 
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corresponds  to  X^.  Equation  (7)  can  also  be  used,  following  the  argument  in 
Ince,^  to  show  that  the  derivative  of  the  characteristic  function,  at  X^,  is 
zero.  Consequently,  X^  is  not  a  simple  zero  of  the  characteristic  function. 
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EXAMPLE 


Let  2i  -  p  e'i  ^3  and  Z2  ■  p  e'  ^  be  a  pair  of  values  given  by  Gros- 
jean  and  De  Meyer ,6  where  p  is  a  positive  root  of 

Ji/3(2/3  p3/2)„0 

and  Ji/3  is  the  Bessel  function  of  order  1/3.  Based  on  equation  (6a), 
the  gradient  is  purely  imaginary  and  is  given  by 


g  .  -8  p3  j 


(10) 


From  equation  (6c),  the  quantity  h  -  X  Is  given  by 


h  -  X  -  4  p3  [  cos(7i/3)  +  sin(7i/3)  i  ] 

H  2 


(11) 


Recall  that  k2(x)  =  g  x  +  h  and  k2(0)  =  h.  Let  h  be  chosen  as 


h  .  8p3  5i=!M3)  i  , 


(12) 


where  c  is  the  nominal  sound  speed  in  m/s  and  f  is  the  frequency  in 
hertz.  The  real  part  of  h  is  the  same  as  the  real  wavenumber  squared. 
The  imaginary  part  of  h  is  designed  to  make  k2(H)  *  {2k  f)2/c2.  The 
multiple  eigenvalue  Xp  is 


X„  ^  .|p3  sin2(7t/3)  cos(7t/3)  ^  ^p3  sin2(7i/3)  ^ 
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The  horizontal  wavenumber  kn  corresponding  to  Xn  is  kn  «  (Xn)^^2> 
where  the  square  root  has  a  positive  imaginary  part.  The  complex  phase 
velocity  is  defined  by  Cp  *  {2%  f)/  kp.  The  eigenvalues  occur  in  a  strip 
(see  the  appendix)  in  the  complex  X  plane  defined  by  |  Im  X  |  ^  a  and 
Re  X  <  b,  where 


a  .  8p3  sin^(-/3) 
and 

b.<2^LQi 

C2 

The  multiple  eigenvalue  in  equation  (13)  is  in  the  middle  of  this  strip. 

The  characteristic  function  will  be  constructed,  and  its  behavior 
at  the  multiple  eigenvalue  will  be  investigated.  A  solution  of  the 
differential  equation  in  equation  (la)  that  satisfies  the  boundary 
condition  u(0,  X)  *  0  in  equation  (1b)  is  given  by 

u(x,X)  -  Ai[z(0.X)]  Ai[e®‘z(x,X)]  -  Ai[z(x,X)]  Ai[e®'z(0,X)]  ,  (15) 

where  6  «  2ji/3  and  z(x,  X)  =  z(x)  is  defined  by  equation  (3).  The  cube 
root  taken  in  equation  (6b)  places  the  arguments  of  the  Airy  functions 
in  the  lower  half  z-plane,  and  the  pair  of  functions  Ai[z]  and  Ai[e®'z] 
are  the  appropriate  linear  independent  set  for  this  region.  The  Airy 
functions  can  be  computed  with  the  algorithm  of  Schulten,  Anderson, 
and  Gordon. 8  The  eigenvalues  are  determined  by  the  roots  of  the 
characteristic  equation  u(H,  X)  «  0  obtained  from  equation  (1c).  The 
eigenvalues  are  the  zeros  of  the  characteristic  function  F(X),  where 


(14a) 

(14b) 
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F(;l)  -  Ai[2(0.x)]  Ai[e®‘z(H.X)]  -  Ai[z(H.>.)]  Ai[e®‘z(O.X)]  •  (16) 


A  specific  example  will  now  be  constructed.  Let  p  be  the  smallest 
positive  root  of  equation  (9)  computed  with^ 

p  =  [  (3/2)  (2.9025862)  ]  -  2.6663527  .  (17) 

Suppose  that  the  sound  speed  is  c  «  1500  m/s,  the  frequency  is  f  =  25 
Hz,  and  the  thickness  of  the  waveguide  is  H  «  200  m.  The  imaginary  part 
of  k2(x)  is  linearly  decreasing  from  a  maximum  at  the  surface  to  zero 
at  the  bottom.  The  maximum  attenuation  a,  computed  with 

a  «  Im  [k(0)]  ,  (18) 

is  6.0896820  dB  per  wavelength,  or  about  0.1  dB  per  meter. 

Figure  1  is  a  plot  of  the  zero-crossing  curves  of  the  real  and 
imaginary  parts  of  the  characteristic  function  F(X)  in  equation  (16). 

The  curves  are  plotted  in  a  rectangular  region  in  the  complex  k2  =  x 
plane.  The  ReF(X)  *  0  curves  are  solid  and  the  ImF(X)  =  0  curves  are 
broken.  The  rectangle  coincides  with  the  strip  defined  by  equation  (14), 
except  that  Re  X  is  bounded  below  by  (27t  f)2/c2.  where  C  =  1700  m/s  is 
a  maximum  phase  velocity.  The  real  part  of  the  phase  velocity  corre¬ 
sponding  to  the  multiple  eigenvalue  in  equation  (13)  is  1542.8188  m/s. 

The  multiple  eigenvalue  is  contained  in  the  rectangle. 

The  zero-crossing  curves  of  the  real  and  imaginary  parts  intersect 
at  a  zero  of  the  characteristic  function.  The  principle  of  the  argument 
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100*Im  (k**2) 


<T> 

cn 


Note:  The  solid  curves  are  the  zeros  of  die  real  part  and  the  broken  curves  are  the  zeros  of  the 
imaginary  part.  The  plot  is  tor  a  region  in  the  complex  *  A  plane  bounded,  in  physical  units,  by  the 
phase  velocities  1 700  and  1 500  m/s  and  attenuation  values  0. 0  and  6. 089682  dB/waveiength.  The 
axes  are  annotated  with  1 00  times  the  real  and  imaginary  parts  of  ^  a  The  characteristic  function 
of  the  Example  has  a  double  zero  in  this  region  of  the  complex  plane. 
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Figure  1 .  The  Zero-Crossing  Curves  Plotted  for  the  Real  and  Imaginary  Parts  of  the 
Characteristic  Function  of  the  Example 


implies  that  two  curves  intersect  at  a  simple  zero  and  four  curves 
intersect  at  a  double  zero.  It  is  apparent  from  figure  1  that  the  first 
eigenvalue,  with  the  smallest  phase  velocity,  is  a  double  zero  of  the 
characteristic  function.  There  will  be  a  generalized  eigenfunction  of 
rank  2  corresponding  to  the  first  eigenvalue.  The  second  eigenvalue  is 
simple  zero  of  the  characteristic  function. 

Evidence  of  the  double  zero  can  also  be  seen  by  considering  a 
sequence  of  values  of  the  parameter  p  that  surrounds  the  value  in 
equation  (17).  These  values  correspond  to  waveguides  with  different 
attenuation  gradients.  The  plots  shown  in  figure  2  are  similar  to  the 
plot  in  figure  1 ,  except  that  p  has  been  given  the  values  2.6,  2.66, 
2.665,  2.668,  2.67,  and  2.7.  In  figure  2(a),  with  p  *  2.6,  there  are  three 
simple  zeros.  In  figure  2(c),  with  p  =  2.665,  the  first  two  zeros  have 
moved  close  to  each  other  and  are  aligned  horizontally.  In  figure  2(d), 
with  p  «  2.668,  the  first  two  zeros  are  still  close,  but  now  they  are 
aligned  vertically.  Finally,  in  figure  2(f),  with  p  *=  2.7,  there  are  again 
three  well-separated  simple  zeros.  The  change  in  alignment  occurs  at 
the  value  of  p  corresponding  to  the  double  zero.  A  movie  made  from  a 
large  number  of  plots  like  those  in  figure  2,  with  a  fine  sampling  of  p, 
would  show  that  the  first  two  zeros  coalesce  and  then  separate  again. 


(b) 


(e) 


(c)  (f) 

Note:  The  parameter  p  has  been  given  the  values  2. 6,  2. 66,  2. 666,  2. 668,  2. 67,  and  2. 7.  The 
sequence  proceeds  from  top  to  bottom  and  left  to  right.  The  double  zero  occurs  between  the  bottom 
of  the  first  column  and  the  top  of  the  second  column.  The  parameter  p  determines  the  attenuation 

gradient  in  the  waveguide. 
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Figure  2.  A  Sequence  of  Plots  for  Values  of  the  Parameter  p  Surrounding  the  Value 

Corresponding  to  the  Double  Zero 


CONCLUSIONS 


The  usual  separation  of  variables  (normal  mode)  solution  for  a  point 
source  in  the  waveguide  of  the  Example  in  the  previous  section  is  not 
valid.  This  is  because  the  usual  normal  mode  solution  assumes  that  the 
poles  of  the  Green's  function  are  simple  and  that  the  residues  at  the  poles 
can  be  computed  by  equation  (8).  The  Green's  function  in  the  Example  has 
a  double  pole,  which  involves  a  more  complicated  formula^  to  compute 
the  residue.  Thus,  the  normal  mode  solution  must  be  modified  to  include 
both  eigenfunctions  and  generalized  eigenfunctions.  A  solution  can  also 
be  obtained  with  an  alternate  technique  like  the  fast  field  program,'^ ^ 
which  computes  the  contour  integral  solution  directly  without  using 
residue  theory. 

It  Is  clear  from  the  derivation  of  the  Example  that  there  are  many 
such  elementary  examples.  The  amount  of  attenuation  required  to  produce 
the  double  pole  in  the  Example  is  quite  substantial.  When  the  amount  of 
attenuation  is  very  limited,  it  has  been  shown  that  all  the  poles  are 
simple."’  ^  It  remains  to  be  determined  if  double  poles  must  be  treated 
systematically  In  normal  mode  computer  codes  operating  in  more 
realistic  undenwater  acoustic  environments. 
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APPENDIX 

BOUNDS  FOR  THE  LOCATION  OF  THE 
COMPLEX  EIGENVALUES 


This  appendix  derives  two  inequalities  that  bound  the  location  of 
the  complex  eigenvalues  of  the  depth-separated  wave  equation. 
Consider  the  differential  equation 


+  [q(x)->-]u  -  0  . 

dx  [p(x)  dx 


(A-1) 


where  q(x)  « r(x)  +  i  s(x)  is  a  complex-valued  function  and  p(x)  is  a 
positive  real-valued  function  on  the  interval  0  <  x  <  H.  It  is  not  initially 
assumed  that  u  satisfies  any  specific  boundary  conditions.  Let  the 
superscript  *  stand  for  complex  conjugation.  If  equation  (A-1 )  is 
multiplied  by  u*  and  the  sum  and  difference  are  formed  with  the 
complex  conjugate  of  the  resulting  equation,  then  it  follows  that 


u*Jl_L_c!y.  ±^iJL  .JL-SLi- 

dx[p{x)dxj  dx  [p(x)  dx. 


[{X±X}  -  (q(x)  ±  q‘(x))]l^ 

p{x) 


(A-2) 


Equation  (A-2)  is  integrated  by  parts  on  the  interval  [0,H]  to  yield 


jiLdu.  ±  -  U-diL- 
p(x)  dx  p(x)  dx  Jo 


1  dul  du  ±  du  dul  dx 
p(x)  L  dx  dx  dx  dxJ 


(A-3) 


-  [  [(A.  ±  A.*)  -  (q(x)  ±  q*(x))]iA^  dx 
Jo  P(x) 


A-2 


If  u  -  <|)n  is  an  eigenfunction  and  X  ■  Xp  is  the  corresponding  eigen¬ 
value,  then  the  boundary  conditions  at  0  and  H  can  be  used  to  eliminate 
the  first  term  on  the  left-hand  side  of  equation  (A-3),  and  it  becomes 


/•” 

r 

1  1 

6^ 

2 

± 

2 

X 

Q. 

0 

- 

dx 

dx 

- 

dx  - 


(A-4) 


where  the  identity  |  z  ■  22*  has  been  used. 


Taking  the  plus  sign,  equation  (A-4)  yields 


I  -i_ 

d<l>n 

L  p(x) 

dx 

dx 


\±l 

P(X) 


dx  -  Re  Xn 


P(x) 


dx 


(A-5) 


Because  the  integral  on  the  left-hand  side  of  equation  (A-5)  is 
positive,  it  follows  that 


Re  Xfi  < 


P(x) 


(A-6) 
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Let  r^ax  “  fTiax{r(x):  0  ^  x  ^  H}.  Replacing  r(x)  in  the  integral  in 
equation  (A-6)  with  its  maximum  gives  the  strict  inequality 


Re  Xn  <  rmax  • 


(A-7) 


Taking  the  minus  sign  in  equation  (A-4)  gives  the  equation 


Let  Sf^jn  -  min{s(x):  0  ^  x  ^  H}  and  Sfy^ax  “  ^ax{s(x):  0  <  x  ^  H}.  Replacing 

s(x)  in  the  integral  in  equation  (A-8)  with  its  minimum  and  maximum 
gives  the  inequality 


Sfnin  ^  Ini  Xn  ^  Smax 


(A-9) 


In  the  depth-separated  wave  equation,  p(x)  «  p(x)  is  density,  which 
is  positive,  and  q(x)  -  k2(x)  is  the  square  of  the  wavenumber.  Thus,  Re 
Xp  is  less  than  the  maximum  of  r(x)  «  Re  k2(x),  which  is  determined  by 
the  minimum  sound  speed.  The  imaginary  part  of  k2(x),  which  causes 
attenuation,  is  non-negative.  Hence,  Im  Xp  is  non-negative  and  less  than 

or  equal  to  the  maximum  of  s(x)  »  Im  k2(x).  This  places  Xp  in  an  infinite 
half-strip  in  the  upper  complex  X  -  k^  plane. 
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